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Abstract. We consider a periodic magnetic Schrodinger operator H , de- 
pending on the semiclassical parameter h > 0, on a noncompact Riemannian 
manifold M such that H 1 (M, R) = endowed with a properly discontinuous 
cocompact isometric action of a discrete group. We assume that there is no 
electric field and that the magnetic field has a periodic set of compact magnetic 
wells. We suppose that the magnetic field vanishes regularly on a hypersur- 
face S. First, we prove upper and lower estimates for the bottom \o(H h ) of 
the spectrum of the operator H h in L 2 (M). Then, assuming the existence of 
non-degenerate miniwells for the reduced spectral problem on S, we prove the 
existence of an arbitrary large number of spectral gaps for the operator H h 
in the region close to \o(H h ), as h — •> 0. In this case, we also obtain upper 
estimates for the eigenvalues of the one-well problem. 



1. Preliminaries and main results 

Let M be a noncompact oriented manifold of dimension n > 2 equipped with 
a properly discontinuous action of a finitely generated, discrete group T such that 
M/r is compact. Suppose that H 1 (A4,M.) = 0, i.e. any closed 1-form on M is 
exact. Let g be a F-invariant Riemannian metric and B a real- valued T-invariant 
closed 2-form on M. Assume that B is exact and choose a real-valued 1-form A on 
M such that dA = B. 

Thus, one has a natural mapping 

u ih du + Am 

from C^°(M) to the space Vl\(M) of smooth, compactly supported one- forms on 
M . The Riemannian metric allows to define scalar products in these spaces and 
consider the adjoint operator 

(ihd + A)* : fl£(Af) — C C °°(M). 

A Schrodinger operator with magnetic potential A is defined by the formula 

H h = (ihd + A)*(ihd + A). 

Here h > is a semiclassical parameter, which is assumed to be small. 

Choose local coordinates X = (Xi, . . . , X n ) on M . Write the 1-form A in the 
local coordinates as 

n 
3 = 1 
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the matrix of the Ricmannian metric g as 

9( X ) = {9jl{X))\<j,l< n 

and its inverse as 

g(X)- 1 = tf'iX))^^. 
Denote |<7(X)| = dct(g(X)). Then the magnetic field B is given by the following 
formula 

dA k dA, 

j<k J 

Moreover, the operator H h has the form 
^9VQ\^<S 9X 



^\gJx)W e {x) (ih^L. + MX)) 



For any x E M, denote by B(x) the anti-symmetric linear operator on the 
tangent space T X M associated with the 2-form B: 

g x {B{x)u, v) = B x (u, v), u, v € T X M. 

Recall that the intensity of the magnetic field is defined as 

Tr+(B(x)) = ]T Xj (x) = Itt([B*(x) ■ B{x)]^). 

\j(x)>a 

i\j(x)€<y(B(x)) 

It turns out that in many problems the function x i— > ft,-Tr + (B(x)) can be considered 
as a magnetic potential, that is, as a magnetic analog of the electric potential V in 
a Schrodinger operator —h 2 A + V. 

We will also use the trace norm of B(x): 

\B{x)\ = [Tr{B*(x) -Bix))] 1 ' 2 . 

It coincides with the norm of B(x) with respect to the Riemannian metric on the 
space of linear operators on T X M induced by the Riemannian metric g on M . 

In this paper we will always assume that the magnetic field has a periodic set of 
compact potential wells. More precisely, put 

6 = min{Tr+(S(a;)) : x £ M} 

and assume that there exist a (connected) fundamental domain T and a constant 
e > such that 

(1.1) Tr + (B(x)) >b + € , xedT. 

For any ei < e , put 

U ei ={xef : Tr + (B(x)) <6 + ei}. 

Thus U €l is an open subset of T such that U ei n dT = and, for ei < eq, U ei is 
compact and included in the interior of T . Any connected component of U ei with 
£1 < e and also any its translation under the action of an element of T can be 
understood as a magnetic well. These magnetic wells are separated by potential 
barriers, which are getting higher and higher when ft — ► (in the semiclassical 
limit). 
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We will consider the magnetic Schrodinger operator H h as an unbounded self- 
adjoint operator in the Hilbert space L 2 (M) and study gaps in the spectrum of 
this operator, which are located below the top of potential barriers, that is, on the 
interval [0, h(bo + eo)]. Here by a gap in the spectrum a(T) of a self-adjoint operator 
T in a Hilbert space we will mean any connected component of the complement of 
<r(T) in R, that is, any maximal interval (a, b) such that (a, b) n cr(T) = 0. 

The problem of existence of gaps in the spectra of second order periodic dif- 
ferential operators has been extensively studied recently. Some related results on 
spectral gaps for periodic magnetic Schrodinger operators can be found for example 
in0[Il[i[Tg[Ti[li[2Ol[2Tl[23[2l[M[23^ (see also the references therein) . 

In our case, the important role is played by the tunneling effect, that is, by the 
possibility for the quantum particle described by the Hamiltonian H h with such an 
energy to pass through a potential barrier. Using the semiclassical analysis of the 
tunneling effect, we showed in [7] that the spectrum of the magnetic Schrodinger 
operator H h on the interval is localized in an exponentially small neighborhood 
of the spectrum of its Dirichlet realization inside the wells. This result reduces 
the investigation of some gaps in the spectrum of the operator H h to the study 
of the eigenvalue distribution for a "one-well" operator and leads us to suggest 
a general scheme of a proof of existence of spectral gaps in [S]. We disregard in 
this paper the analysis of the spectrum in the above mentioned exponentially small 
neighborhoods. 

We consider the case when bo = and the zero set of the magnetic field has 
regular hypersurface pieces. More precisely, suppose that there exists xo £ M such 
that B(xq) = and in a neighborhood U ofxo the zero set of B is a smooth oriented 
hypersurface S, and, moreover, there are constants k £ Z, k > 0, and C > such 
that, for all x £ U, we have: 

(1.2) C^dix, S) k < \B(x)\ < Cd(x, S) k . 

On compact manifolds, such a model was introduced for the first time by Mont- 
gomery [26] and was further studied in [TUl EH El [9] . 

Denote by N the external unit normal vector to S and by N an arbitrary ex- 
tension of N to a smooth vector field on U. Let Wo,i be the smooth one form on S 
defined, for any vector field V on S, by the formula 

(V,u , 1 )(y) = ^N k (B(N,V))(y), y £ S, 

where V is a C°° extension of V to U. By (|1.2p . it is easy to see that wq x{x) ^ 
for any x £ S. Denote 

Wmin(-B) = inf \w ,i(x)\ > 0. 

For any a £ R, denote by Ao(a) the bottom of the self-adjoint second order 
differential operator 

d 2 ( t k+1 \ 2 

in I/ 2 (R). Put v := inf Qe R Ao(a). It is clear that v > 0. We refer the reader to 
Section [2] for more properties. 

In [9], we have proved the following result. 
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Theorem 1.1. For any a and b such that 

v w min (B) < a < b 

and for any natural N, there exists ho > such that, for any h G (0, ho], the 
spectrum of H h in the interval 

2k + 2 2k + 2 

[h' T ^ T a, h~ k ~^ r b] 

has at least N gaps. 

In this paper we will concentrate our analysis on the region close to the bottom 
of the spectrum. First, we state upper and lower estimates for the bottom Xo(H h ) 
of the spectrum of the operator H h in L 2 (M) 

Theorem 1.2. Suppose that the operator H h satisfies condition (|l.ip with some 
€q > 0, and that the zero set of the magnetic field B is a smooth oriented hypersur- 
face S. Moreover, assume that there are some k 6 Z, k > and C > such that, 
for all x in a neighborhood U of S, we have: 

C- l d{x, S) k < \B(x)\ < Cd(x, S) k . 

Then there exist C > and ho > such that, for any h G (0, ho), we have 

, ■ 2 2k + 2 6k + 8 , , 2 2k + 2 6fc + 8 

i> u>m in (B) k +z h k + 2 - Ch^ k + 2 ~> < X (H ) < v uj min (B) k + 2 h fc + 2 + Ch 3 ( k + 2 > . 

A similar result was obtained for the bottom of the spectrum of the Neumann 
realization of the operator H h in a bounded domain in E 2 by Pan and Kwek [55] 
in the case k — 1 and by Aramaki [T] in the case k arbitrary odd. 

As an immediate consequence of Theorems 11.11 and 11.21 we obtain the following 
statement. 

Corollary 1.1. In addition to the assumptions of Theorem \1. 6 A suppose that M 
is compact. Denote by X (H h ) < Xi(H h ) < X2(H h ) < ... the eigenvalues of the 
operator H h in L 2 (M). Then, for integer m > 0, we have 

lim h~^X m (H h ) = v w min (S)^T5 . 

h^0 

In the case when k = 1 and |wo,i(a;)| is constant along S, this result was obtained 
by Montgomery [26) . In Theorem 1.4, under some additional assumptions, more 
generic than in |26j , we obtain stronger estimates for the eigenvalues of the one- well 
problem. 

Like in the case of the Schrodinger operator with electric potential (see [16]). 
one can introduce an internal notion of magnetic well for a fixed hypersurface S 
in the zero set of the magnetic field B. Such magnetic wells can be naturally 
called magnetic miniwells. They are defined by means of the function |wo,i| on S. 
Assuming that there exists a non-degenerate miniwell on S, we prove the existence 

, 2fc+3 , 

of gaps in the spectrum oiH on intervals of size h k + 2 , close to the bottom Xo(H ). 

Theorem 1.3. Suppose that the operator H h satisfies condition (jl.ip with some 
eo > 0, and that there exists Xq G M such that H(xq) = and in a neighborhood 
U of xq the zero set o/B is a smooth oriented hypersurface S. Suppose moreover 
that there are constants k G Z, k > and C > such that for all x G U we have: 

C^dix, S) k < \B(x)\ < Cd(x, S) k . 
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Assume finally that there exists x\ G S such that |cJo,i(a;i)| = Wmin(-B) and, for all 
x G S in some neighborhood of X\ 

Cid s {x,xi) 2 < \w ,i(x)\ -u> mia (B) < C 2 d s {x,xi) 2 . 

Then, for any natural N , there exist > and h^ > such that the spectrum of 
H h in the interval 

2 2 fc -|- 2 2 2fc-|-2 2 fc -|- 3 ~ 

i> w m in(B) *+2 fe"^ , i> w min (B) w /iW +b N h~*+z 
has at least N gaps for any h G (0, /ijv)- 

As an immediate consequence of Theorem 11.31 we also obtain upper bounds for 
the eigenvalues of the one- well problem. 

Theorem 1.4. Suppose that M is a compact manifold of dimension n>2. Assume 
that the zero set of B is a smooth oriented hypersurface S, and there are constants 
k G Z, k > 0, and C > such that for all x G S we have: 

CrHip, S) k < \B(x)\ < Cd(x, S) k . 

Assume that there exists x± G S such that |wo,l(^i)| = ^min(-B) and, for all x G S 
in some neighborhood of x\, we have the estimate 

Cid s {x,x 1 ) 2 < |wo,i(x)| -uj min (B) < C 2 d s (x,x 1 ) 2 . 

Denote by Xo(H h ) < \\{H h ) < \2(H h ) < ... the eigenvalues of the operator 
H h . Then, for any natural m, there exist C m > and h m > such that for any 
h G (0, h m ) we have 

\n{H H ) < i>W min (B)^2/lW# + C m h^r. 

Under the additional assumption that there exists a unique miniwell, we believe 
that, using the methods of [4j, one can prove the lower bound for the ground state 
energy X (H h ) of the form 

, 2 2fc + 2 2k + 3 

Ao(# ) > j>w min (B)^T2/i~ _ Ch~^ 
and the upper bound for the splitting between Xo(H h ) and \\{H h ) of the form 

\i(H h ) - \ {H h ) < Ch^. 
For this, we have to know a certain property of the ground state energy \o(a) of 



the operator (|1.3j) . which we state as a conjecture. 

Conjecture 1.1. Any minimum of Ao(oi) is non-degenerate, that, is, for any 
G K such that Ao(o; m in) — v we have 



(1.4) ^^ mi ") > °- 

Moreover, we believe that one can prove the lower bound for the splitting between 
X {H h ) and A^iJ' 1 ) of the form 

X 1 (H h )-X (H h )>Ch^, 
if, in addition, the following conjecture is true. 

Conjecture 1.2. There exists a unique a m i n G K such that Ao(a m in) = v. 
We refer the reader to Section [2] for discussions. 
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2. Some ordinary differential operators 

For any a £ K and (3 £ M., f3 ^ 0, consider the self-adjoint second order differen- 
tial operator in L 2 (R) given by 

* ' 1 r #w- a ' 2 



V 



In the context of magnetic bottles, this family of operators (for k = 1) first appears 
in [25] (see also [TO])- Denote by \q(cx, (3) the bottom of the spectrum of the operator 
Q(a,/3). 

In this section, we recall some properties of Ao(a,/3), which were established in 
[2"B1 HU1 [25] , First of all, let us remark that Ao(a,/3) is a continuous function of 
a E R and (3 g M. \ {0}. One can see by scaling that, for (3 > 0, 

(2.1) A (a,/3) = f3*h\ {[3~^a,l) . 

A further discussion depends on the parity of k. 

When k is odd, Ao(a, 1) tends to +00 as a — > — 00 by monotonicity. For analyzing 
its behavior as a — * +00, it is suitable to do a dilation t — aws, which leads to 
the analysis of 



a 



ds 2 



-.k+l 



1 



with h — a ( fc + 2 )/( fc + 1 ) small. One can use the semi-classical analysis (see [3] for 
the one-dimensional case and [301 HI] f° r the multidimensional case) to show that 

A (a, 1) ~ (k + l)^n"a*rr j as a +00 . 

In particular, we see that Ao(a, 1) tends to +00. 

When k is even, we have Ao(a, 1) = Ao(— a, 1), and, therefore, it is sufficient to 
consider the case a > 0. As a — > +00, semi-classical analysis again shows that 
Ao(a, 1) tends to +00. 

So in both cases, it is clear that the continuous function Ao(a, 1) is nonnegative 
and that there exists (at least one) a m i n € M such that Ao(a, 1) is minimal: 

A (a min , 1) = i>. 

The results of numerical computation^ for a m i n , v and the second eigenvalue Ai 
of the operator Q(a m i n , 1) are given (modulo 10 -2 ) in Table [TJ 



k 


1 


2 


3 


4 


5 


6 


7 


^min 


0.35 





0.16 





0.10 





0.07 


V 


0.57 


0.66 


0.68 


0.76 


0.81 


0.87 


0.92 


Ai 


1.98 


2.50 


2.61 


2.98 


3.18 


3.47 


3.66 



Table 1. Numerical results for a m i n , v and Ai 



^performed for us by V. Bonnaillie-Noel 
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In Fig. [T]and [2j one can also see the graphs of the function A = Xo(a, 1) and its 
quadratic approximation at a = a m i n : 

A quad (a) = A (a min , 1) + l|^(a min , l)(a - a min ) 2 . 



Figure 1. k even 



Figure 2. k odd 

Numerical computations show that when k is even the minimum is attained at 
zero: a m i n = 0. They also suggest that the minimum a m j n is non degenerate, 
supporting Coniecture ll.li and that the second derivative ^jr(a m i n , 1) tends as k 
tends to oo to 2. One more step towards the proof of Conjecture [T7T] can be done 
by the following considerations. 

Let u° £ I/ 2 (R), || = 1, be the L 2 normalized strictly positive eigenvector of 
the operator Q(a, 1), corresponding to the eigenvalue \a(a, 1): 

J2 o / f k+i \ 2 
(2.2) Q(a, 1)< = -^f + (— a) < = A Q (a, 1)<. 

One can show that u° a depends smoothly on a. Differentiating (|2.2j) with respect 
to a, we obtain 

(2.3, Q (a , - 2 (£L - „) < = ^(„, 1,< + ««, 
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that implies that 
(2.4) 
and 
(2.5) 



da 



(a,l) 



d 2 \ 



da 2 

It follows from (l2~4l) that 
(2.6) 

and, for k odd, 



(a,l) =2-4 



t k+i 
k + l 



— -a] «(t)) 2 dt 

k+l a() J da 
(u° am Jt)fdt = 0, 



fk+l 

m «,„W) 2 ^>o. 

It has been claimed that this minimum is unique for k = 1 in [28] and for arbitrary 
odd fc in p. 

It also follows from (|2.3[) that 

fk+l 

(2.7) ' ' ' "" ^ ' ' 



du° 

(Q(a min , l)-u)^ = 2 
da 



k + l 



Note that the above computations can be made not only for a minimum point 
a m j n , but also for any critical point of Ao. 

We will also need the following identity (see [2"8] , Proposition 3.5 and the formula 
(3.14)): 

2 



(2.8) 



1 



k + l 



_ t k+i 



a 



min I "'a., 



k + 2 



Based on numerical computations, one can give a proof of Conjecture 11.11 for 

dvP _____ 

small k as follows. Since is orthogonal in L 2 (R, dt) to u° a , by (|2.7p and (|2.8[) . 



we have 



< 



1 



< 



Ai — v 
2 

Ai — v 
2 



l)-i>) 



9a ' 9a 



fe+i 



fc + 1 



This implies that 



Ai - v \k + 2 / 
dui 



^1/2 


dul 




da 



da 



\i-v \k + 2 



1/2 



and 



/ 



2v 



iTT-°)<-^"* s (5_--»)(t + j) 



By (|2~5]) . it follows that 



S// 



da2 [a min , i) _> (Ai _j> )(fc + 2) 



= 2 



(fc + 2)A_ - (fc + 6)£> 
(fc + 2)(A_ -i>) 
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Hence if for some k, we have 

(2.9) 0+2)Ai > (fc + 6)i>, 

we deduce that the corresponding minimum a m i n is non-degenerate. 

Using the data given in Table [H one easily verifies that the condition (|2.9[) is 
satisfied for k — 1, . . . , 7. Due to the accuracy of these numerical computations, 
one can consider that (| 1 .4|) is safely controlled for k = 1, . . . , 7. 

When k is odd, the potential of the Sturm-Liouville operator Q(a, 1) is even and, 
by well-known properties of Sturm-Liouville operators, any eigenfunction associated 
with the m-th eigenvalue A m (a) of Q(a, 1) is even if m is even and odd if m is odd. 

In particular, vP a is even for any a, and, therefore, -^f is even as well. This shows 
that one can replace Ai by A2 in the above arguments, stating that if, for some odd 
k 7 we have 

(fc + 2)A 2 > (fc + 6)z>, 
the corresponding minimum a m i n is non-degenerate. 



3. ESTIMATES FOR THE BOTTOM 

In this section, we will prove Theorem 11.21 So we assume that the zero set of 
the magnetic field B is a smooth oriented hypersurface S, and the magnetic field 
satisfies the estimate (|1.2I) in a neighborhood U of S with some k E Z, k > 0. 

Let G be the Ricmannian metric on S induced by g. Denote by dxQ the corre- 
sponding Ricmannian volume form on S. Let 



i* s A 



be the closed one form on S induced by A, where i$ is the embedding of S to M. 
For any t G R, let Pg ^Wo.o + ^i^ k+ltJ o,ij be a formally self-adjoint operator in 
L 2 (S,dxa) defined by 



x ( ihd + wo.o 



1 



t ^0,1 



Consider the self-adjoint operator H h '° in L 2 (R x S,dtdxa) defined by the for- 
mula 

d 2 „, ( 1 



H 



h .0 



= -It 



dt 2 



+ Pg> I W 0j0 + 



k + 1 



fe+i 



^0,1 



By Theorem 2.7 of [10], the operator H h '° has discrete spectrum. 
The quadratic form associated with the operator H h '° is given by 



q h [u] 



r+00 p 


h 2 


du 


2 

+ 


/-oo Js 




dt 





(ihd + ujq.o + 



1 



9o 



dt dx„ 



u e C C °°(IR x S). 



Without loss of generality, we can assume that U coincides with an open tubular 
neighborhood of S and choose a diffeomorphism 



9 : I x S -» U, 
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where I is an open interval (—so, £o) with sq > small enough, such that |{o}xs = 
id and 

(3.1) (e* 3 -g )| {0 }xs =0, 
where go is a Riemannian metric on I x S given by 

go = dt 2 + G. 

By adding to A the exact one form d<j>, where <p is the function satisfying 
N(x)(p(x) = -{N,A)(x), xeU, 
4>(x) = 0, i£S, 

we may assume that 

(3.2) (N,A}(x) =0, xeU. 

Denote by Hp the unbounded self-adjoint operator in the Hilbert space L 2 (D) 
defined by the operator H h in the domain D = U with the Dirichlct boundary 
conditions. The operator Hp is generated by the quadratic form 

u !->• q D [u] := / \(ih d + A)u\ 2 dx 
Jd 

with the domain 

Dom(g^) = {ue L 2 (D) : {ihd + A)u G L 2 Q}{D),u\ eD = 0}, 

where L 2 Vl}(D) denotes the Hilbert space of L 2 differential 1-forms on D and dx is 
the Riemannian volume form on D. Denote by Xo(H D ) the bottom of the spectrum 
of the operator H D , By Theorem 2.1 in [7], there exist C,c,ho > such that for 
any h G (0, ho] we have 

\\ (H h ) - \ (H*)\ <Ce- c '^. 

It can be seen from the proof of Theorem 2.7 in [TU] that, if there exist h > 0, a 
family of functions w h G C%°(I x S), h G (0, ho], and a function A°(/i) defined on 
h G (0, /io] such that 

x (h)<dh (2k+2 '>^ k + 2 \ he(o,h ], 

and 

||(fr h ' - \°{h))w h \\ < C 1 hP k+a M h+a '>\\w h \\, he (OM], 

with some positive constants C and C\, then there exist h\ < ho and a positive 
constant C*2 such that 

||(#£ - X"(h))v h \\ < C 2 h (2fe+3)/(fc+2) ||v' l || , Vh G]0, /ii] , 

where u' 1 = (0~ 1 )*u>' 1 G C^°(U) . So, in order to complete the proof, it sufficient 
to prove that there exist h > and positive constants C\ and C*2 such that, for all 
h &]0,h ], 

x 2 2fc + 2 Sfc + 8 , n 2 2fc + 2 6fc + 8 

^Wmin^)^ 2 ^ 1 ^ -Ci/l 3 < fc + 2 ) < A (-ff < zJWmh^-B)'^ 2 /! 1 ^ + C 2 /l 3 < fe+2 > . 

The desired upper bound for Xo(H h, °) follows from the construction of approximate 
eigenfunctions of H h '° given in [5]. It remains to prove the lower bound for \o(H h '°). 
We will localize the problem in two scales. 
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First, we choose a covering of S by local coordinate charts, S = \J 



m—l 



u„ 



Let Xm € C°°(S) be a partition of unity subordinate to this covering so that 
supp Xm C U m for each m and 



X 2 m (x) = 1, xeS. 



m—l 

There exists a C°° real valued function <p m such that on U„ 

^0,0 = d(p m . 
Therefore, for v e C£°(U m ), we obtain 

q h [v] 



+ 



rs 


h 2 


dv 


l-oo JS 







(ihd+ -7——t + w ,i) exp --<p. 

fc + 1 \ Al 



Secondly, for a fixed m = 1, 2, . . . , d, by scaling a standard partition of unity 



pn— 1 



we can find a partition of unity satisfying 



Ei4wi 2 = 1 . 



and 



suppy 



C 



»,3 



:= {x e 



\X Vm,j 



j\<£oh P }, 



where (3 > is a parameter which will be determined later. 

Choose an arbitrary z m ,j £ Qm,j- Let g m j — G(z m ,j) be a flat Euclidean metric 



pn— 1 



For any I and p, we have 



|G" p (x) - gXj(x)\ = 0(hP) as h - 0, a: e Q m ,,, 



(3.3) - 
We also have a similar estimate for Wo,i: 

(3.4) \uo,i{x) - uo,i(z m j)\ = 0(h f3 ) as h -> 0, ieQ 



Consider the self-adjoint operator ° in L 2 (K x Q m j , dt dx 9mj ) defined by the 
formula 



H 



h,0 



<9t 2 



The quadratic form associated with H^f- is given by 



/+oo ,■ 
\ 
-oo J Q 



1 


h 2 


dw 






dt 



+ 



{ihd+ Y^t k+1 GJ 0tl (z mJ ))w 



\G(z mj )\^ 2 dtdx, 



w G C C °°(1R x Q m>j ). 
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Lemma 3.1. For any w G C^°(R x Q m ,j), we have 

qtjw] > f>\uj 0tl {z m ^)\T^hW\\ w \\ 2 . 
Proof. Using the simple scaling t = h}/^ k+2 \, we can write 



•> Qm,j 



dw 



dr 



+ 



1 



(^ fc + 2 d+ j-j-jT k+1 u; 0A (z m j))w 



\G(z m!j )\^ 2 drdx. 



There exists a (? m> j-orthonormal frame {X\, X 2 , ■ ■ ■ , X n -i} in R n 1 such that 

Wo,i(z m> j){Xi) = \u> ,i{z m ,j)\, uJa,i(z md )(Xi) = 0, I > 2. 
Then, by a linear change of variables in R n , we obtain 

dw 1 







dw 


2 


ri 

1 —OO J Qm,j 




dr 


+ 



2k + 2 

>h 



1=2 
r-DC 



dw 



n 



dr alyx dy 2 ■ ■ ■ dyn-! 



dw 



dr 



■ , -J— 9w 1 u , 1 . . . . 



Denote by w(j\\,y 2 , • • • , JM-i) the partial Fourier transform of w in the yi-variable. 
We have 



/ 





dw 


2 




dr 


+ 



dw 1 



lhk + 2 + __^+l|^ 1 ( Xmj )| w 



drdy 1 dy 2 ...dy n - 1 



dw 



dr 



k + l 



|wo,i(2 m> j)| - h k + 2 rji) \w\ drdr}idy 2 . ■ -dy n -i. 



For any fixed r]i, y 2 , ■ ■ ■ , y n -i, the expression 



dw 



dr 



+ 



1 



-T k+i \uJo.i{z m . J )\-h^r] 1 j \w\ 2 



dr 



is the quadratic form of the operator Q(a,f3) with 

a = h^r) 1 , P = \uJo,i{z m j) 
evaluated on id. Therefore, we have 



dw 



dr 



1 



-T k+1 \uo,i(z m ,j)\ - /j- fc + 2 m J \w\ 



k + 



dr 



/OO ^OO 
\w\ 2 dr > v\aJo,i(z m ^)\T^2 / |w| 2 dr, 
-OO J —OO 



that immediately completes the proof. 



□ 
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On the other hand, using the inequality 

2\ab\ < e 2 a 2 + e- 2 b 2 , 

with e = hP (here p > is a parameter which will be determined later) and estimates 
(pT3]l and IpO]) . we get, for any w £ C^°(R x Q m ,j), 



~>]=(i + o(^)) 



/■+OO / 




ft 2 


9w 


/ — oo </ £ 






0* 



1 



T^-y* w ,i)(exp ( ~-<f m ) w 



\G(z mJ )\ 1/2 dtdx 



>(l-^)(l + 0(^))g^.[exp(-- % 



/i- 2p (l + 0(^)) 



1*1 



2fc+2 



x |(cj 0l i(a;) - u) ,i(z m ,j))w(t,x)\ 2 . |G(z mj -)| 1/2 df cfe 



>(1 - h 2 ")(l + 0(//)X 3 [exp f--^ 
-_^*-*(l + 0(^)) 



|i| 2fe+2 |w(t,a;)| 2 |G(a;)| 1/2 Adx. 
Therefore, using the IMS formula, for any v E C£°(R x Z7 m ), we obtain 

^m=E^[^]-^Eiii v ^hi 2 

> [(1 - h 2 P) (1 + 0(^))i) Wmill (B)w^ - ch 2 - 2 ? 
and, furthermore, for any u <E C^°(]R x S), we get 

d d 



q h [u} = 52q h {Xmu}-h 2 Y^ IIIVX, 

m— 1 m= 1 



> [(1 - ft 2 ") (1 + 0(//))Pu, min (S)^ft 



(fc + 1) ; 



:h 20-2p( 1 + 0(hP))\\\t\ k+1 u\ 



It is shown in [10] (see Remark 5.2) that, given some constant C, then there 
exists, for any k e N, a constant such that, for ft G (0, 1] and any eigenvalue 
X°(h) of if' 1,0 such that 

A (ft) < a/i (2fe+2)/(fe+2) , 
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then any corresponding eigenfunction U h satisfies 

|||t| fe+ V|| <C k hV' +1 M k+ Q\\u h \\. 
Therefore, we obtain, for some new constants Co, Ci, C2, C, 

X°(h) > (*>aw(B)^ - Coh? - Cih 2 ? - C 2 /i 2 "- 2 ') hW - Ch 2 - 2 ? . 
Putting (3 — 3 (k +2 ) > P = 3(fe+2) ' ^his thieves the proof of the lower bound. 

4. Existence of gaps 
In this section we will prove Theorem II .31 



4.1. A model operator. First, we will only assume that there exists xq E M 
such that 'B(xq) — and in a neighborhood U of xq the zero set of B is a smooth 
oriented hypersurface S, and the magnetic field satisfies in U the estimate (II. 2[) 
with some fe 6 Z, ft > 0. As above, we will assume that U coincides with an open 
tubular neighborhood of S and choose a diffeomorphism O : I x S — > U, where / is 
an open interval (— £o,£o) with eo > small enough, such that |{o}xS — id- We 
will also assume that the conditions (|3.ip and (|3.2p hold. 

Suppose that 5 admits a coordinate system with coordinates (si, . . . , s n _ 1). 
Then we have a coordinate system in J7 with coordinates X = (Xq, X\, . . . , X n _i) 
with Xq — t G I and X,- = Sj, j = 1, 2, . . . , n — 1. Thus, S 1 is given by the equation 
t = 0. We have 

Vm-X-)] <a,/3<n-l 

where 



v£=i/i^ + A a pO, a = 0,l,...,n-l. 
or, equivalently, 

(4.1) H h = V g«P{X)V h a V h + ih V r«v£, 



0<Q,/3<n-l 0<a<n-l 



where 

r° = 1 - V ( J\gTx)\g^ a (X)) , a = 0,1, 

By (|3~T|) and (f3T2]) . we have 



A (t,s) = 0, 

and 

. 9 = (it 2 +G + 0(f), <^0, 
where G is the induced metric on S. So we can write 
g 00 (t,s) = l+g m {s)t + 0{t 2 ), 
go 3 {t,s)^g 0j (s)t + 0(t 2 ), 
gjl (t, s) = G jt (s) + g 3 i(s)t + 0(t 2 ). 

In particular, we have 

\g(t,s)\ = \G(s)\+0(t), t^O. 
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By assumption, we can write 

Mt, a) = ^fa) + ^4%)^ + ^- 2 ^l( s )t^ + 0(t k +% t - . 

Suppose that a family u h £ C^°(U),h 6 (0, /io], satisfies the following assump- 
tions. 

For any real m > 0, there exists C m > and h m > such that, for any h 6 (0, /i m ), 
we have 



(4.2) 

(4.3) 
and 



\t m V h 3 u h \ 



\\t m u h \\ <C m h^\\u h \\, 

h\\t m ^\\<C m h^\\u h 
at 



(4.4) \\t m \7^u h \\ + h ( ||i m ^V^|| + \\t m V'}—u h \\ 



d 2 u h 



h 2 \\? 



m + 2fc + 2 ,, . ,, 

< C m h-i+*-\\u h \\ 



As shown in |10j . a family u h 6 C%°(U), h S (0, /io), such that, for some C > 0, 
we have 

(H h u h ,u h ) <Ch^T\\u h \\ 2 , he(0,h o ), 

satisfies the conditions P~2")) , (TO]) and P~4")l . 
By (|4~Tj) . it follows that 



— V" + V" — 



Of 



dt 



du h 



dt 

Observe that the first and the third terms can only contribute to the terms of order 

2fc + 2 

h k + 2 . More precisely, we have 



H h u h 



2„,h 



,d 2 u 
dt 2 



l<j,£<n-l 
X 



Y g^( S ) (ih± + + kh^)t k+ 1 



This fact was stated in \W\ Theorem 2.7] . For the proof of Theorem ll.3l we have 

2k + 3 

to improve the remainder 0{h k + 2 ) in the last formula. For this purpose we will 

h h • 1 

take into account further terms in the expansion of H u in powers of h k + 2 . This 
leads us to introduce a new model operator Hq given by 



HI 



& 



d 2 



dt 2 



dt 2 



(8)1* 



fc+1 



d_ 

dt 



l<j<n-l 
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+ £ G^( S )^A + W 0) (s)+ i « (s)t ^i 

l<j,i<n-l 



fc + 2 



l<j/<n-l V S J 

x(*^+^w + ^4!!(»)*«) 

l<?<n— 1 ^ ^ ! 



U \s)t k+1 



Lemma 4.1. Suppose that a family u h € C%°(U), h G (0, ho), satisfies the condi- 
tions (|4.2[) . (|4. 3[) and (14. 4[) . Then, there exists C > suc/i i/ia£, /or any /i € (0, /io), 
we /iaue 

HffV* -H^u h \\ < Ch 2 \\u h \\ . 

4.2. Approximate eigenfunctions: main result. Now we additionally assume 
that there exists x± £ S such that |wo,i(xi)| = w m in(B), a neighborhood V of x\ in 
S and a constant C\ > such that, for all ieV, 

Cid s (x,xi) 2 < |w ,i(x)| -^min(-B) < Cid s {x,xi) 2 . 

Take normal coordinate system / : U(xi) C S — > K™ -1 on 5 defined in a neighbor- 
hood E/(xi) of xi, where f(U(xi)) = B(0, r) is a ball in R" -1 centered at the origin 
and / (xi) = 0. As above, we will denote local coordinates by s — (s\, s 2 , . . . , s n _i). 
Note that 

1/2 

^)=lEi^(o)p' 

Consider the Euclidean space K™ _1 with coordinates (o"i, 02, • • • , o"„-i). Take 
the unit vector 

^ = — ^ (45(o), ■ ■ ■ ,4?r 1} (o)) e k"- 1 , 

W m in (.-£>) V y 

and complete it to an orthonormal base (ei = e w , ei, . . . , e„_i ) in K™" 1 . Denote by 
e^,, e2, . . . , e n —i the corresponding first order differential operators with constant 
coefficients in In particular, we have 

The Laplacian A = — ~§^i m K™ -1 can be written as 

A = A w + A w i , 

where 

n-1 

A w = -e 2 , A w x=-^ef. 

i=2 
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Consider a second order differential operator K in W 1 1 given by 

i d 2 \ , 



(4.5) 
where 



K 



2 da 2 



(amin, f )A W + A w i + fi rm oyer m + A , 



2) 5^: (0) + a ^ 1- ^7 (0) ^r (0) 



and 



(0)a„ 



tl — 1 p h-\-1 k-\-l 

+ 2c min (B)- 2 £4>)4>) / l_(I_- amin )« min ( T ))2 dr 
+ ^ min (B)- 2 ]T .9^(0)4^(0)^(0) 

l<j,£<n-l 

x /r(^-a min ) « ml „(r)) 2 dr. 

It should be noted that we have indeed the operator K = K ainin attached at any 
minimum a m ; n and we can do the same construction at any a m i n - 

A construction of approximate eigenfunctions of the operator Hp in D = U is 
given in the next theorem. 

Theorem 4.1. For any critical point a m j n and for any A in the spectrum of K = 
K ainin , there exist Cy > 0, hi > 0, and a family U h £ C£°(D), h £ (0,/ii], such 
that, for any h £ (0, hi], we have 



{H h D -z(h))U h \\ L ,<Cih^\\U h \\ L 2, 



with 



z(h) = i>uj nlin (B) k + 2 h k + 2 + Xh fc + 2 . 

Proof. The proof of this theorem is long and will be divided into several steps. □ 

4.3. Formal expansions near the minimum. Choose a function <p £ C°°(B(Q, r)) 
such that d<f> = cjq,o- F° r some a £ 1, we make a gauge transformation 



u(t, s) = exp 



Then we have 



As) 



exp 



( J j }(0) Sj \ 

■ 3=1 

1 fc + i l 

Wmin(B)^/l fc + 2 



v(t,s) 



seB(0,r), t£R. 



(4.6) H£ u(t, s) = exp 



exp 



■ 3 = 1 

1 fc+i l 



P h v(t,s), 
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where 

p h = p h + p h + p h + p h + p h + p h i 

and 

P}=2ih «"(»)' 

l<j<n-l 

P 3 fc =ift £ 5° J M^(s)i fe+1 , 

l<j<n-l 

Pf = E G ^ S )( Z/l ^- + fcTT^ (S)tfe+1 ~ «"nun(B)-S^42(0) 
l<j,£<n-l J 

- aW m in (P) ~ Wg^i (0) + -i-^^t^ 2 ), 



V ^ (s)t ( jft ^ + ^_„ g(8)tW , 

l<j,£<n-l J 

- au, min (P)-^/i^«(0)) (iftA + ^-^(5)^+1 

- au min (By — h—uj^l (0) J , 

p 6 h =-h 2 n(s)§- t -h e r o( s )(^+^r^](^ fe+i 

l<j<n-l J 

fc+l fc+l / '\ \ 

- aw min (B)-*nfc*T3w$(0)J . 
We now make the change of variables 

(4.7) t = ^ min (B)-^2 &V(fe+2) T) g = h^+^a , 

and expand the operators in powers of h as h — > up to the terms of order 

Q l {Ak+7)/{2(k+2))y p or the firgt three termS) 

we obtain 

2k + 2 , . 2 d 2 2k + 3 , , 1 nn , , 9 2 , 4fc + 7 

P 2 h =2ihW ^'(0V min (P)-^o;g(0)r - a) |- + 0(*«), 

l<j<7l-l ^ / T 

l<j<n-l 

For the analysis of the fourth term, let us start with the computation of 

' =( ih Wj + kh J " { ' )tl *' - ™™»(sr™'' ia <4:!(°) 
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d 



1 



fc + 1 fc+1 



i 



k + 2 °- 2V ' 
After the change of variables, we obtain 



^ 2fc+2 2fc + 2 / T 

J=W mln (S) <= + 2 /l-S+2" 



fe+1 



^CO^/a^a)^ _ C0 (0) 



x(^ TT 4!l(^ 1/2(fc+2) -)-a4!ko)) 



+ 



2fc + 3 



d I r 

3 

d 2 



r k+l 



d<7jd<7£ 



2fc + 3 T 



9 
<9t7£ 



/c + 2 



x ( (irT^ (/ll/2(fe+2)<7) ~ aw S (0) j 4^ 1/2(fe+2 M 
+ (^4> 1/2(fc+2 M -a45(o)) 4^ 1/2(fe+2 M 



We can write 



(4.8) u%l(hW**>v) = cg(0) + ftW&T £ ^1(0)^ 



<9s r 



9 



2, 



2 ds r ds m 
Since s = is a minimum of the function 



(0)a r a m + (3(/i 2 «=+ 2 > ), /i— >0. 



(4.9) 



we have 



(4.10) 

Using (|4.8p . we obtain that 



T k+l 

X>)t 



^(ftVWV) _ ^) (0) j f 1^45 (^M ff) _ aw W (0) 
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2 



=45(0)43(0) - a) 2 + ^ (41(0) £ ^ ( 0)«x 



fc+i 



-fe+i 



fc + 1 \k+l 



d 2 J l) 



E 



— a 



r fc+i 



1 V fe + 1 



<9s T 



-(0). 



(0)<T r CT„ 



r 2fc+2 



(fc+1) 2 



Similarly, we have 



_3 

<9(7j V & + 1 



43(/> 1/2(fc+2 M-*43(o) 



fc+1 



(9 



+2^^^ K 



+ /j2(fc + 2) 



9cty 



fc + 1 



+ 0(/l^F2). 



Thus, we get the following expression for /: 



f=o; min (B)-^/ l ^4fl(0)43(0) (J 

2fc+2 i 1 
+ /l"E+2"" l "2(fe + 2) 



— a 



fc+1 

(s) -ffi h o, (0) ^ +4! , (0) ^_ 

*4l 



fc + 1 



+ c, min (i?)-W#(^)(0)E9j i (0K 

+4S(o)E^(°w 

a 2 



fc+l V fc+ 1 



2fc + 3 
+ /l <= + 2 



dojdoi 



+ iu m[n (B) *=+2 



(0) 



r fc+i 



do^ fc+1 



fc+2 . _fc+l 



- min (i?)-w# (45(0)43(0) +43(0)4^(0)) raGbrH 



+- min (s)-^[i(45(o)E^f-( K-. 



3 2 4 j) 



+4>)E 



(0)aycr TO ) 



-fc+i 



r fc+i 



fc+1 Vfc + 1 



— a 
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1^^— (O)-aT- 



9s r <9s m (fc + l) 2 

Using (|4.10|) and the expansion 



4fc + 4 
+ 0(/l^+2j). 



d G 3 3 

(0)cr r cr m + 0(/l 5 ^+ 3 T), fc->Q, 



2 ds r ds n „ 

we obtain that, after the change of variables (|4.7[) . the fourth term takes the form: 

r fc+i 



Pi 1 =u; min (B)^h^- 



2fc+2 / T 
+ 2 



fc + l 

+ ft,"E+2- + 3^T2) 2io; m in(S) _ '^t 



fc + l 



E<l(o) 



2fc + 3 



d 2 



EU , „. fc+l \ - 



(9(7, 



(0) 



fc + l 



r fc+2 . _fc+l 



2o; r 



EE 



~2k+2 



fe+1 / -fc+l 



fc+l Vfc+ 1 



<9s r 



(O)^i(Q) a r a n 



ds 



(fc + l) 2 



ds r ds T 



-(0)45(0)45(0) ^ 







-l \ 2 - 














Vfc4 







+ 5£ E 

4fc + 4 

Finally, after the change of variables (|4.7p , the fifth and sixth terms become: 



^^(fl)"W E ^(^)4fl(0)4!l(0)r 

l<j,£<n-l 



fc+l 



4fc + 4 

+ 0(/i 2 (<=+ 2 >), 



P 6 h = -^(0)^(5)^^^- 

- E ^^ # r J (OV min ( J B)-^4fl(0)(^— _a) + 0(h^T). 
l<j<n-l 

Thus, after the change of variables (|4.7p . the operator P^ has a formal asymptotic 
expansion 

oo 

■—1 , > 2 2fc + 2 ^ — > 1 ^ 

P = W min (P)^+2-/l"^+2- \ h^+^Pi, 



22 



B. HELFFER. AND Y. A. KORDYUKOV 



where 

d 2 ( r k+1 

(4.11) P = __+f___ a ) =Q(a,l), 



(4.13) 



(4-12) Pi = 2i Wmin (B)-^ _ a j J2<m d(y / 

and 



-k+l 







2 



d_ 

dr 2 

2* £ ^(0Mni„(S)-^^](0)r(^-a)| r 



+ * ^ ^"(0)<](0)^ m i„(S)-^r fc+1 

l<j<n-l 



— a 



2a, mm (B)-«S ! Y>g(0)4»(0)^(l^ 

» («)- 2 [E(E»S(o)|i^(»))-...^(^ 

r,m j v 



a 



<9s r ds m ) (k + 1) 2 



3 

2 



a 



2 

^(0)4fi(0)4>)r ( ^ - a 1 

l<j,£<ra-l X " 

T°(0)a; min (S)-^ — 



Wmin (s)-* 53 rj(0Kl(0) 

l<j<n-l 



jfe + 1 



4.4. Reduction to the zero set. Now we use the method initiated by Grushin 
[5] (and references therein) and Sjostrand 31J in the context of hypoellipticity. We 
will closely follow the exposition in [4] (see also [6]). We now choose some a m i n 
and will use the previous construction at 



Oi Q^min ■ 

pn— l-i 



The starting point is to consider the operator in <S(R n ) x 5(M™ ) defined by 



Pa Ro 
R+ 
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V = Q(a m ln, 1) - V, 



where the operator Pq : <S(R") — » 5(R n ) is given by 
d 2 ( r k+1 

the operator : S(W 1 ^ 1 ) -> S(R n ) is given by 

i?o 4>{t, a) = </>(cr)< ml „ (r), 4> G 
and the operator : S(R n ) -> ^(R"" 1 ) is given by 

Jtf/M = / fir, aK<T)dT, f e <S(R"). 



We observe that Vq considered as an operator in L 2 (R ra , dr da) x L 2 (R™ _1 , da) is 
formally self-adjoint. In particular Rq is the Hilbertian adjoint of Rq (considered 
as an operator from L 2 (R™, dr da) to L 2 (M n_1 , da)). 
We also verify that 

where n : L 2 (R") — > i 2 (R") is the orthogonal projection on the subspace {Rw° . }x 
L 2 (R™ _1 ) in L 2 (R«): 

n /(r,a) = H /(r,,)< mm (T)drj < mlI >), / £ L 2 (R"). 

Define the operator Eq in 5(R") by 

(4.14) e = (i -n Q )p^\i -n Q ), 

where, by abuse of notation, we consider Pq as an operator in L 2 (R"). As shown 
in [4j Lemma A. 5], Eq respects the Schwartz space ^(R™). Then we have 

VqoSq = I, 

where the operator £q in 5(R™) x <S(R ra_1 ) is given by the matrix 



So 



Eq R q 

R+ 



The idea is to consider the more general operator V(z) in 5(K") x S(M. n 1 ) 
defined by 

V{z) = ( ^ain(B)-^2h- 2 -^T(P h -Z) R^ 



l 

R£ 



Note that Viz) is for z G R formally self-adjoint for the original L 2 -scalar product 
in i 2 (R") but not for the usual L 2 associated to the standard Lebesgue measure 
dr da). 

We are looking for the right inverse of Viz) for any small z G C. One can write 

m + » So" 

where 

(4.15) SP = uj min (B)-^T2h-^TTP h - Pq 

= -0 + h'^+^Pi + h~^P 2 + h^A^Q(h), 

(4.16) Z = w miu (B)~T^h~^z - v, 
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and Q(h) admits a complete expansion 



We first observe that 
where 

K 



3=0 

V(z)£ = I + K, 

(SP - Z)E a (SP - Z)R+ 




We will assume that Z is a function of h, Z = Z(h), which admits a formal 
asymptotic expansion of the form 

(4.17) Z(/i)~^2,/iW. 
Then we have 

SP- Z ~ ~ z t) h '^ TT) ■ 

l>\ 

If we define 

q ~ E(- 1 ) J/C ' > 

3=0 

then the operator is well-defined (after reordering) as a formal expansion in powers 
of h^js+aj anc j 

V(z)£ Q Q~I. 
So £(z) = £qQ is the right inverse of V(z). If we write 

E(z) E+(zf 
.E-(z) E ± (z) / 



£{z) 



we get, in the sense of formal expansions in powers of /i 2 ( fc + 2 > , 



(4.18) ^{B)-^h-^-{P h - z)E(z) + RqE-{z) ~ / 



2fc + 2 



(4.19) w mln (S)-^2/i- ^tr{P h - z)E+(z) + R^E ± (z) ~ 

(4.20) RqE(z) ~ 

(4.21) R+E+(z)~I 

Let us introduce a function E ± (Z) by 

E±(Z) = E ± (z), 
with Z related to z by (|4.16p . We have 

Ki ^ ([{hp - z)E Q y [(sp - z)E Q y- 1 [(sp - z)r 



and therefore, E ± (Z) is the following asymptotic series, 

(4.22) E±{Z) ~ ]T(-1)^ + [(£P - Z)E y- 1 [(SP Z)R- 



3 = 1 
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Lemma 4.2. The function E (Z(h)) admits the following formal asymptotic ex- 

i 

pansion in powers of /i 2 < fc + 2 ) : 



with 



^{Zih)) = Y E± 3 h ^ 



i=i 



(4.23) 
(4.24) 



E i — Zi, 

2(fc + 2) 



E\ = Z 2 -u min (B)- — K . 



Proof. It follows from (I4.22p that the coefficient of ft, 2 ( fc + 2 > is given by 



2(k + 2) 



Using the definitions of the operators Rq , P\ and R Q and (|2.6[) . we get 



i - 2ia; min ( J B)- — 

2(fc + 2) 



r k+l 



k + 1 



"min ) |Uq .(t)| 2 o!t 



By (|4.22p . the coefficient of h k + 2 is given by the operator 



(4.25) 



E i — Z 2 — R^P2R + Rq PiEqPiRq 



By (j443]> . it follows that 



R£P 2 R^ = UJ min (By— ^ q^2 + X! ^m^r^m + ' 

j 3 r,m 



where b rm and a are given by 



rm — ^min 



Rn 



k + 1 \k + l 



Q2 / ,,(j) 



«.-(E^!(0)^(0)) 



(k + 1) : 



2 "« U+l 



fc+i * - 



«.-|E|5-(»vS(0)4'!(0) 
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and 
(4.26) 

a = - g m (0)u; min (B)-—> / r^f (r)< mIn (r)dr 



k+ „ r / T k+1 \ a„.0 

+ 2« E 5° J (0V m in(B)- — ^(O)yr^— -a min J-^(rK min (r)dr 



l<i<n-l 



A r fe+1 

fe2(^TT" amin ) (u °- (r))2dT 



/■ / fe+1 \ 

+ c min (S)-W# £ ^( S )u;g(0)4>) /r(^— -J « m ,>)) 2 dr 



i f fin 

-r»c min (£0-— y ^(r)< mm (r)d 



r 

-u; min (i?r^l ^ rj(0) W gl(0) f f^--n irAlt )i,r jr.,,^/. 

By (123)) and ((231) , it follows that 
+ ( rk+1 \ rk+1 



/ T fe+1 \ r fc+l 

( vm -a mi „J m « mm (r)) 2 dr = 



fe + 2' 



2 

min ' 



T 2fc+2 r T 2k+2 £ 

and 



^O" T—r - "min i? = II _ "min U L, n ( r )l| 2 = 



fe + i y u 11 Vfc + i / min fc + 2 

Using the fact that Gjg(0) = Sji and dGji/ds r (0) = for any j, £ and r, the 
formula (|4.9[) implies that, for any r and m 

f^(°)=E^(°)4>)4>) 

as r ds m *rt os r ds m 

^E^<oh<;!(o) +2 E^(»)^(»)- 

From the above formulae, it follows that 

brm = ~ kWn(-E0 fc + 2 ^rm- 

Using integration by parts and (|2.6| . we get 
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that implies that the sum of the second and the third terms in (|4.26[) equals zero. 
It is also easy to see that the last two terms in (|4.26p equal zero. Thus, we have 

(4.27) a = -u> min (B)-T& A. 
We conclude that 

(4.28) R^P 2 Ro 

r,m 

Next, by flTi"]) and (|4~T2]) . we have 
E PiR 

j_ _ / r fe+i \ 

= 2iuj min (B) {I-U )(Q(a min , l)-z>) 1 {I-H ) ( T^jy - "mm J < mln (r)e w . 

By ([23]) . it follows that 



n 



(fTT ~ amin ) u °™ (r) = / (jTi " amin ) (u °™ (r))2dr = °- 



Therefore, by (|2~7)) . we get 

r fe+i 

- Ho) ' ' 



/ T fc+i \ / r fe+i \ 



du° 

(Q(a mln ,l)-i>)^. 



^0 

We obtain (note that ilo-^f = 0) 



i du° 

E PiR^ = iw„in(B)-srs-=-aeL 



Next, we have 

k + 1 / oa 



du° 

R+P 1 E P 1 Ro = 2tu min (B)-^R+ [ — — - a mln ) — ^A u 



T k+l 



By ([2~5]) . it follows that 

T k+1 \ Oil" 



Rt 



\k + l mll V da 

= /(* + !- H U ^ (T) ^ (T) " T = 4 ( 2 - ftF^' 1} 
Therefore, we obtain 

(4.29) R+P 1 E P l Ro = iu; min (£)-W5 ^ - ^(a min , 1)^ A„. 

From (|4~25]) . (|4~28|) and ([4~29]) . we get ![3~24"|) . that completes the proof of the 
lemma. □ 
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4.5. Construction of approximate eigenfunctions. In this section, we com- 
plete the proof of Theorem 14.11 So suppose that A is in the spectrum of the 
operator K. Our considerations depend on whether Conjecture [TTT] is true or false. 
Since this is unknown at the moment, we consider both possible cases. 

First, suppose that Coniecture ll.ll is true, that is, the second derivative ^r(a m in, 1) 
is positive. Then the operator K has discrete spectrum. Let 4>o 5(M n_1 ) be an 
eigenfunction of K with the corresponding eigenvalue A. Put in (|4.17|) 

Z t =0, Z 2 = u mia (B)-*T2\, Z £ = 0(W>3). 

So we have 

2 1 . 2 2fc + 2 1 2fc + 3 

Z (h) = uj min (B) k + 2 Xh k + 2 , z(h) = i)uj m i n (B) k + 2 h k + 2 + Xh k + 2 K^ 2 . 
By Lemma 321 we have 

E^iZQi))^ = 0{h^) ) , 

and, by (14.19|) . we obtain 

(P h - z{h))E + {Z{h))<j> a = O(h^), 
where the function E + (Z) is given by 

E+(Z) = E+(z), 

with Z related to z by (|4. 16[) . From (|4.6|) . it follows that the function 



U h (t, s) = X (t, s)h-^J exp 



exp 



( a ffim E41(0)^ 



V / 

xE + {Z{h))<j) {ij min {B)^h- x ^ k+2 \h- l ^ k+2 h\ s€B(0,r), t e K, 
where x G C£°(f ) i s a cut-off function, satisfies \\Uh\\ = 1 + o(l) and 

(H' l - z(h))U h = Oih^rh). 
From the above, we see that E + (Z) is the following asymptotic series, 

E+(Z) ~ i? - + ^2(-iyE [(6P - Z)E ] j - 1 [(SP - Z)B^\, 

i 

and, therefore, it admits an asymptotic expansion in powers of h 2 <- k + 2 > : 

oo 

E + {Z(h)) ~V£^_/iW5). 

— ' 2(fc + 2) 
l=\ 

Using this fact, one can easily see that the function U h satisfies the conditions 
P~3]) and (|44|) . By Lemma |4~TI it follows that 



{H h - z{K))U h = 0{h^+v), 
that completes the proof in this case. 
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Now consider the case when Conjecture II. II is false. So suppose that the equality 
^r(oirain, 1) = is true. Thus, the operator K has the form 

K = A L0 ± + 22 ^rm^rO'rn + A. 



Lemma 4.3. There exists Wq £ S(W l 1 ). \\wq\\ = 1, which satisfies the following 
conditions: there exists C > such that, for any h > 0, we have 

\\(K-\)w%\\ < CTi 5 ^, 

and for any multi-index a — (oti, . . . , a n —i), there exists a constant C a > such 
that, for any h > ; we have 

(4.30) ||fl£u# || < C a h~^+V. 

Proof. Take an arbitrary vector e' u £ K n_1 , which is orthogonal to the vectors 
e%, ■ ■ • , e„_i with respect to the positive definite bilinear form fl in R n_1 given by 

fi(<7, </) = ^2 ^rm,0'rO J m- 

Consider the linear coordinate system in R n_1 with coordinates {pi,p2, ■ ■ ■ , p n -i) 
defined by the base (e^, e<2, ■ ■ ■ , e„_i). In these coordinates, the operator K has the 
form 

n-l q 2 ™~ 2 
j=2 "l r,m=2 

where Vt' n > and the quadratic form X)™ m 2 =2 ^'rmPrPm is positive definite. 

With respect to the decomposition L 2 (M™ _1 ) = L 2 (R P1 ) ® £ 2 (R™~ 2 Pn _ t )), the 
operator K can be written &sK = V®I + I® K$, where V is the multiplication 
operator in L 2 (M P1 ) by the function 

V(jh) = n' llP i, Pl £ R, 

and Kq is the harmonic oscillator in L 2 (R p l ~ 2 Pn _ 1 ) given by 

n— 1 n— 2 

^0 = - ^ 7f2 + Y ^'rmPrPm- 

Therefore (see, for instance, [29, Theorem VIII.33]), the spectrum of K equals 



a(K) = a(V) + a(K ) = [*o(K ), +oo), 

where Xq(Kq) — inf ct(Kq). 

Denote by uq £ iS(K"~ 2 ), ||uo|| = 1, the cigenfunction of Kq associated with 
Ao(A'o). Define the function w'q £ S(M n ~~ 1 ) by the formula 

(4.31) w h Q {p) = ch- 1 '^e-^-^l^ hmk+2) \o{p2, Pn-i), P E R-\ 

where the constant c > is chosen in such a way that \\w'q\\ = 1 and a satisfies the 
condition 

A = X Q (K Q ) + f^'na 2 ■ 

Then 
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+ C fc-Va(^) n / i ^ e -0. 1 -a)V(2^^) Uo(pa) . . . ) pn _ l)) p e R n-1 ; 

and, therefore, we have 

< Ch^ k+2 l 

The estimates (|4.30|) follow immediately from the explicit formula (|4.31|) for Wq . □ 
Take Wq as in Lemma POl and 



The operator — (Pq — v) 1 (I — IIo)^p2 — Pi-EoPlJ has the form of a second order 

differential operator in cr, whose coefficients are bounded operators in L 2 (M., dr). Us- 
ing this fact, it can be easily checked that, for any multi-index a = (ai, . . . , a„_i), 
there exists a constant C a > such that, for any h > 0, we have 



(4.32) 
Put 



<9>? < C a h~ 



v h (T,a) =(i?o -h^2jE P 1 R^)w^{T,a) + h^2R^w'l(T,a) 

i i Bit 

oa 

+ h^W^T, cr). 

2 2fc-|-2 2 h | 3 

Then, with #(/i) = z>w min (i?)"fc+2 + \h~^ , we have 

-"v. , 2 2fc + 2 r 1. 

(P' 1 - =u; min (P)^^ [(P - 

+ fr^W (Pxi?-^ - (Po - t>)P PiPo «>o) 

+ ^ (((P 2 - w min (P)-^2 A)P.,7 - PiPoPiPq )^o 

+ (P -i>)wf +0(/i 5 ^ 5T ) 



2 2fc + 3 
=W min (P)* + 2 /l fe + 2 



n ((P 2 -ow(P)-^A)Po 



--h—R^(K - A)wq + 0{h^^), 



Therefore, we obtain 



||(P /l - z{h))v h \\ = 0(h^+^). 
From (14.61). it follows that the function 



[J h (t, s) = x{t, s)h 4 " fe + 2 > exp ( -i 



exp 



V 



n-1 



ttmin E W J l(°) S i 

fc+i l 
Wmin(P) fc+2 /l fc + 2 



/ 



x /(WminiB)^^ 1 ^^ 2 ^, /T 1 / 2 ^ 2 ^), s G P(0, r), t € K, 



PERIODIC SCHRODINGER OPERATORS WITH HYPERSURFACE MAGNETIC WELLS 31 



where y is a cut-off function, satisfies \\Uh\\ = 1 + o(l) and 

(H$ - z(h))U h = 0(h$&). 

Using (|4.30[) and (|4.32[) . one can easily verify that the function U h satisfies the 
conditions P~2")) . (|Q1) and By Lemma [JU it follows that 

(H h - z(h))U h = O(h^k) , 

that completes the proof of Theorem 14.11 

If Coniecture ll.ll is true, following the arguments of @J, one can prove the follow- 
ing refined version of Theorem 14. II (cf. 4, Theorem 3.1]): for any A in the spectrum 
of K, there exist a sequence {Q}^_ ® ano - a sequence of functions {0j}j*L o m 
C%°(D) such that, for any TV > 0, there exists M > such that, if 

M 

. . , 2 2fc + 2 2fc + 3 4fc + 7 x — , j 

z M (/i) = vbj^(B)T^h~^ T + Xh~^ + /i^fe+2j \ hw^jQ . 



and 



M 

3=0 

then 

|| (H h - z M (ft)) ^ < 0(h N )\\cj> h M \\L*, h - 0. 
Moreover, if, in addition, Conjecture 11.21 is true and a miniwell x\ G S is unique, 
one can show, following the lines of [4], that the m-th eigenvalue X m (H h ) of the 
operator H h admits the asymptotic expansion 

oo 

, 2 2fc + 2 2fc + 3 4fc + 7 ^ — ^ j 

X m (H ) = vuo min {B)T^h-^ +X m h— + h^"+^ 2^/i 2 < fc+2 >0, 

3=0 

where X m is m-th eigenvalue of the operator K. 

4.6. Proof of Theorem II .31 To complete the proof of Theorem ll.3l we will use a 
general result on the existence of gaps in the spectrum of the magnetic Schrodinger 
operator H h on the interval [0, h(bo + eo)] obtained in [8, Theorem 2.1]. Fix t\ > 
and 62 > such that e\ < £2 < £0j an d consider the Dirichlet realization on 
the operator H h in the domain D = U €2 . The operator Hjj has discrete spectrum. 

Theorem 4.2. Let N > 1. Suppose that there is a subset fig < ^\ < .. . < ^ of 
an interval 1(h) C [0, h(bo + e±)) such that 

(1) There exist constants c > and M > 1 such that 

l x^-^_ 1 >ch M , j = l,...,N, 

dist(4, 81(h)) > ch M , dist(/4, 81(h)) > ch M , 

for any h > small enough; 

(2) Each Hj,j = 0, 1, . . . , N, is an approximate eigenvalue of the operator H^: 
for some v'j G C^°(D) we have 

\\Hbv$-P$v$y i =a j (h)\\v$\\, 
where ctj(h) — o(h M ) as h — > 0. 
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Then the spectrum of H h on the interval 1(h) has at least N gaps for any sufficiently 
small h > 0. 

Fix any N > 1. If ^jr-(a m in, 1) > 0, then K has discrete spectrum Ao < Ai < 
. . ., Xj — > oo as j — > oo. Take an arbitrary 6jy > Ajv- By Theorem 14.11 for any 
j = 0, 1, . . . ,N, there exist C) > 0, h 0J > 0, <pj(h) G Cf(D) and Zj(h) with 

2 2k + 2 2fc + 3 4fc + 7 , 

Zj(/i) = i>u! m i n (B) k + 2 h k + 2 + Xjh k + 2 + 0(h 2k ~+z) 
such that, for any h G (0, /lo.jL we have 

|| (Hb-z i (h))<t>(h)\\<C J h%&\\<l>(h)l 
So we apply Theorem 14.21 with 

f 2 2fc + 2 2 2fe + 2 2fc + 3 

1(h) = 0uj min (B) k + 2 h k + 2 , i> uj min (B) k + 2 h k + 2 + b N h~ k + 2 

and /ij = Zj(h),j = 0, 1, . . . , N, that completes the proof in this case. 

If ^jT-(a m i n ,l) = 0, then the spectrum of if is a semi-axis [Ao,oo). Taking 
an arbitrary on > Ao and arbitrary Ai < . . . < Xn on the interval (Ao,6/v) and 
proceeding as above, we complete the proof. 
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